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In piezoelectric ceramics the changes in the phase stabilities versus stress and tem- 
perature in the vicinity of the phase boundary play a central role. The present 
study was dedicated to the classical piezoelectric, lead-zirconate-titanate (PZT) ce- 
ramic with composition Pb(Zr .54Ti . 46)03 at the Zr-rich side of the morphotropic 
phase boundary at which both intrinsic and extrinsic contributions to piezoelectric- 
ity are significant. The pressure- induced changes in this two-phase (rhombohedral 
R3c+ mono clinic Cm at room temperature and R3c+P4mm above 1 GPa pressures) 
system were studied by high-pressure neutron powder diffraction technique. The ex- 
periments show that applying pressure favors the R3c phase, whereas the Cm phase 
transforms continuously to the P4mm, which is favored at elevated temperatures due 
to the competing entropy term. The Cm — » R3c phase transformation is discontinu- 
ous. The transformation contributes to the extrinsic piezoelectricity. An important 
contribution to the intrinsic piezoelectricity was revealed: a large displacement of 
the B cations (Zr and Ti) with respect to the oxygen anions is induced by pressure. 
Above 600 K a phase transition to a cubic phase took place. Balance between the 
competing terms dictates the curvature of the phase boundary. After high-pressure 
experiments the amount of rhombohedral phase was larger than initially, suggesting 
that on the Zr-rich side of the phase boundary the monoclinic phase is metastable. 
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I. INTRODUCTION 



Piezoelectric lead-zirconate-titanate [Pb(Zr z Tii_ x )03, PZT] solid solution system was 
developed over 40 years ago yet attempts to understand its properties continue to trigger 
new studies. A long-lasting view is that when x is approximately 0.52, a first-order phase 
transition occurs between tetragonal and rhombohedral phases, resulting in two-phase co- 
existence. The electromechanical properties peak slightly on the rhombohedral side of the 
phase boundary. In the composition-temperature plane the boundary (commonly called 
as the morphotropic phase boundary, MPB) is nearly independent of temperature, thus 
making PZT very practical material for applications^. The commonly offered reasoning 
for the exceptionally good electromechanical coupling is based on the idea that there are 
eight (rhombohedral phase) and six (tetragonal phase) spontaneous polarization directions 
available in the two-phase system so that the system can readily respond to external electric 
field or stress. 

The space group symmetries given for a disordered solid-solution should be taken as av- 
erage symmetries from which short-range order deviates. For instance, it has been known 
for long that Raman scattering data cannot be explained by the average symmetries. The 
high-temperature cubic phase has no first-order Raman modes yet experiments revealed 
that spectra collected on PZT above the Curie temperature have rather strong features at 
energies close to the low-temperature first-order phonon energies. In the case of so-called 
relaxor ferroelectrics this type of behavior is normal and the frequently offered explanation 
is that symmetry-lowering defects generate polar nanoregions (see, e.g., refs. 2-^). Also 
the low-temperature Raman spectra of Ti-rich PZT have many features which are not con- 
sistent with the tetragonal symmetry: the twofold degenerate ^-symmetry modes of the 
tetragonal PZT were split, indicating that the symmetry is lower than P4mm£. Raman 
experiments showed that anharmonicity plays a significant role in lead titanate, the anhar- 
monic contribution being increased with increasing temperature 6 . The traditional view was 
modified once high-resolution x-ray synchrotron studies revealed that the phase believed to 
be tetragonal possesses monoclinic distortion^ in the vicinity of the MPB. Neutron powder 
diffraction experiments, able to resolve the monoclinic splits, ruled out octahedral tilts, and 
verified the Cm symmetry^. 

Accurate modeling of the system requires not only the consideration of the unit cell but 



also crystallographic twins (or ferroelectric domains) and grain boundaries must be taken 
into account. In piezoelectric ceramics the response to external stress or electric field can be 
divided into intrinsic and extrinsic contributions 10 . The former is essentially a single crystal 
response (i.e. is formed by the ion displacements within a primitive cell of the crystal), 
whereas the latter covers the contribution due to grain boundaries, preferred orientation or 
texture of the grains, i.e. ferroelectric domains within the grains, and changes in crystal 
phase fractions. Since the full model considering contributions from atomic scale up to 
the macroscopic grain size scale is very complex, experimental studies have commonly been 
applied to gain deeper insight. 

Non-180 domain switching (i.e., contributing to the extrinsic contribution) gives rise to 
approximately 34% of the measured coefficient of PZT 11 . The extrinsic contribution can 
be larger or smaller if the domain wall motion is respectively made easier or more difficult by 
dopingi 2 -^. A study of the domain switching showed that the 90° domains in single phase 
tetragonal phase (titanium rich PZT) hardly switch, whereas the domains in the two-phase 
region switch^. Texture and strain analysis of the ferroelastic behavior of Pb(Zr .4gTio.5i)03 
by in-situ neutron diffraction technique showed that the rhombohedral phase plays a sig- 
nificant role in the macroscopic electromechanical behavior of this material 15 . The domain 
nucleation and domain wall propagation are central factors limiting the speed of ferroelectric 
polarization switching^ 1 ^. 

An important intrinsic contribution to the piezoelectricity is due to the increase of certain 
piezoelectric constants once the phase transition is approached. This increase was predicted 
to be significant in the vicinity of the pressure-induced phase transition in lead titanate^. 
The computations carried out for lead titanate further show that it is the competition 
between two factors which determines the morphotropic phase boundary^. The first is 
the oxygen octahedral tilting, favoring the rhombohedral R3c phase, and the second is the 
entropy, which in the vicinity of the morphotropic phase boundary favors the tetragonal 
phase above 130 K. If the two factors are in balance over a large temperature range, a steep 
phase boundary results in the pressure-temperature plane which is desirable for applications. 
The advantageous feature of the R3c phase is its ability to be compressed efficiently by tilting 
the oxygen octahedra, in contrast to symmetries prohibiting oxygen octahedral tilting 20 . 

We briefly summarize the relationship between the structural parameters and polyhedral 



tilts and volumes, given in ref.—. We follow ref. 



23l and parametrize the asymmetric unit of 



FIG. 1. The R3c phase, whose hexagonal unit cell is shown in panel (a), and the Cm phase, 
panel (b), behave very differently under applied pressure. The Va/Vb ratio between the oxygen 
octahedral and cuboctahedra-volumes of the R3c phase decreases with increasing pressure: the 
crystal is contracting and thus the B cations (which fit oxygen octahedra tightly) have to take 
larger relative volume from the total volume (from the cuboctahedra, which has excess of space for 
Pb) by tilting oxygen octahedra. The symmetry prohibits this mechanism in the P4mm and Cm 
phases. Density-functional theory computations predict that PAmm has an entropy term benefit 
at elevated temperatures. Two rhombohedral (corresponding to the R3m phase) pseudocubic cells 
are shown by dotted lines in panel (a). Due to the octahedral tilting, indicated by arrows, the 
two cells are not equivalent: the tilting corresponds to the Rim — > R3c symmetry lowering. The 
primitive cell of the Cm phase is shown by dotted lines in panel (b). Structure figure was prepared 
by the VESTA software^. 

the R3c phase as given in Table HI 

There is one short and one long 0-0 octahedral edge length parallel to the hexagonal 
afr-plane, labeled as I — Al and I + AZ, respectively (see also Fig. [TJ. Now, the octahedral 
tilt angle is given by tan a; = 3 1//2 4e and the polyhedral volume ratio Va/Vb is equal to 
6 K 2 cos 2 oj — 1, where K is given by equation a = 2Klcosu. The present study focuses 
on the two-phase, Cm and R3c, PZT ceramic material, Pb(Zr .54Tio.46)0 3 , which has a 



TABLE I. The asymmetric unit of the R3c phase as defined in ref. [23. 



x 



y 



z 



Pb 












Ti/Zr 







t 








composition slightly on the Zr-rich side of the morphotropic phase boundary. The main 
goal was to determine the phase fractions and structural parameters as a function of applied 
pressure and temperature. Also the question concerning the reversibility of the structural 
properties of PZT is addressed. 

II. EXPERIMENTAL 

To address the possible homogeneity differences due to the variation in solid-state reac- 
tion based sample preparation method lead zirconate-titanate powders were prepared using 
different starting oxides and sintering conditions. In the first route the PbO, ZrC>2 and 
TiC>2 oxides were mechanically mixed in desired proportions, whereas in the second method 
PbTiC>3 and PbZrC>3 powders were used as starting chemicals. The phase purity and crys- 
tal structure were checked by X-ray powder diffraction and scanning electron microscopy 
measurements. No significant differences were observed and thus a sample prepared through 
the latter method was used for the experiments. Samples were annealed by first forming 
perovskite structure at 1073 K (30 minutes), then increasing the temperature to 1373 K 
(60 minutes) to improve the sample homogeneity and then cooling the sample first to a 
stepwise manner to room temperature. Annealing times were kept rather short in order to 
limit PbO loss. High-pressure neutron powder diffraction experiments were carried out at 
the Los Alamos Neutron Scattering Center using the TAP-98 toroidal anvil press^ 1 ^ set 
on the high-pressure-preferred orientation (HIPPO) diffractometer— >22. Pressure was gen- 
erated using the high-pressure anvil cells. Sodium chloride was used as a pressure calibrant 
material. To minimize deviatoric stress built up during room-temperature compression on 
the polycrystalline sample, all data in our high P-T neutron-diffraction experiment were 
collected during the cooling cycle from 800 K at each desired loading pressure. Data were 
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collected between 300 and 800 K as a function of pressure. Rietveld refinements were carried 
out using the program General Structure Analysis System (GSAS) 28 and EXPGUI 29 . The 
pressure was estimated from the reflection positions of the NaCl phase through the equation 
of stated At higher pressures it was necessary to include the reflections from the diamond 
anvils in the refinement model. The broad hump seen in the background intensity between 
2 and 3 A is due to the diffuse scattering from the amorphous zirconium phophate gasket 
and was modelled using the diffuse scattering option available in the GSAS software. 



III. RESULTS AND DISCUSSION 

a. Structural model. The X-ray diffraction pattern collected on Pb(Zro.54Tio.46)03 pow- 
der is characteristic to the morphotropic phase boundary composition, the most apparent 



indication of a two-phase co-existence is seen 
the R3c + Cm structural model (see refs. [31 



rom the pseudo-cubic 200-reflections. Thus, 



331 ) was used for the refinements of the low- 



pressure data at ambient temperatures. Refinements indicated that the monoclinic distortion 
continuosly vanished with increasing hydrostatic pressure and increasing temperature. The 
monoclinic structure became tetragonal and was correspondingly modelled by the P4mm 
space group. Fig. |2] shows the pattern collected at 3 GPa pressure at room temperature 
and the computed intensity. At ambient conditions the majority phase was monoclinic, see 
the GPa datum in Fig. [3j With increasing pressure the situation changed significantly 
(Fig. [3]), accompanied by drastic changes in rhombohedral tilts and polyhedral volumes (Fig. 
H]). Slight increase of the tetragonal phase fraction with increasing temperature at constant 
pressure is seen in Fig. [3j The lattice parameters given in Fig. H] indicate that the Cm 
phase does not continuously transform to the rhombohedral phase: the difference between 
the rhombohedral and monoclinic structures remains large up to the point at which the 
Cm phase continuously transforms to P4mm phase. Instead, through the studied pressure 
and temperature range yet there are significant changes in the phase fractions. This is in 
line with the first-order phase transition and shows that no continuous polarization rotation 
occurs. Thus, the phase stabilities as a function of pressure and temperature follow well 
the predictions based on the first-principles studies carried out for PbTiOsr^ 1 ^. Further, 
the entropy term seems to have a crucial role for setting the boundary between the pseudo- 
tetragonal and rhombohedral phases: the pseudo-tetragonal phase fraction increases with 
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FIG. 2. Observed (red) and calculated (green) time-of- flight neutron powder diffraction data and its 
difference curve between measured and computed curves (purple) for a Pb(Zro.54Tio.46)03 sample 
at 303 K and 3 GPa. The tick marks, from down to up, are from the R3c, Cm, NaCl (pressure 
standard) and graphite (from the pressure chamber) phases. The statistical figures of merit were: 
X 2 = 2.300, R wp = 2.04 %, R p = 1.42 % and the background substracted R parameters were 
R bwp = 2.75 % and R bp = 1.65 %. 



increasing temperature at constant pressure. 

b. Octahedral tilting. Figure |5] shows the octahedral tilts in the R3c phase and the two 
characteristic octahedral edge lengths, I — Al and I + Al. The octahedral tilt increases with 
increasing pressure, though the tilt angle saturates at high pressures. Thus with increasing 
pressure the volume fraction of the octahedra increases, consistently with the idea that, 
when compared to the tightly filled oxygen octahedra, lead ions have excessive space inside 
cuboctahedra formed from 12 oxygen atoms. In addition to the oxygen octahedral tilting 
also another mechanism can be seen: the continuous expansion of the I + Al and contraction 
of the I — Al. Fig. [6] (a) shows the location (Zr or Ti) and oxygen bond lengths in the 
rhombohedral phase. At ambient conditions the B cations are closer to the larger oxygen 
triangle, consistently with the earlier data^i. This situation changes with increasing pressure: 
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FIG. 3. Rhombohedral weight fraction at ambient conditions and as a function of temperature at 
approximately 1 and 3 GPa pressures. 
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FIG. 4. Lattice parameters of the R3c, Cm, PAmm and Pm3m phases at ambient conditions and 
as a function of temperature at approximately 1 and 3 GPa pressures. Monoclinic and tetragonal 
6m axis lengths are surrounded by a square. The CM-axis values are enclosed by a circle. The Cm 
phase transformed to the PAmm phase at around 400 K at 1 GPa pressure. At ambient conditions 
the monoclinic angle /3 was 90.01(97)° and at 1.26 GPa pressure /3 was 90.62(4)°. The 3 GPa data 
is indicated by a dotted line. Due to the thermal pressure, the pressure values of the highest two 
temperatures (cubic phase) are larger. The inset shows the rhombohedral angle a. 



it is seen that the 5-cations are closer to the small oxygen triangle, indicating that at higher 
pressures the S-cations favour to form a small tetrahedron rather than being centered closer 
to the octahedron center, see the inset of Fig. |6j Positions in which the B cations are closer 
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FIG. 5. Octahedral tilt angles (a), octahedral edge lengths (b) and polyhedral volume fractions of 
the R3c phase at ambient conditions and as a function of temperature at approximately 1 and 3 
GPa pressures. 



to the large triangle is clearly unfavourable as it would result in bond lengths failing to fullfill 
the bond-valence criteria. At 3 GPa pressure the distance between the vertex of the large 
oxygen triangle and triangle center alone is slightly larger than the given B-0 lengths. For 
piezoelectricity this has important consequences: if stress is sufficiently strong, it switches 

10 



the position of the B cations from a larger oxygen triangle towards the smaller oxygen 
triangle thus contributing to the intrinsic piezoelectricity. Thin film technology allows a 
deposition of selected crystal planes in which the biaxial stress can be adjusted by choosing 
the substrate and composition so that the piezoelectric proeprties can be optimized. 

Fig. [6] (b) gives the distance between the oxygen triangles, D(l + Al,l — Al). Figs. Mb) 
and [6] (b) show that whereas D(l + Al, I — Al) and I — Al both decrease and I + Al increases 
significantly when pressure increases from to 1 GPa, D(l + Al, l — Al) hardly changes when 
pressure increases from 1 GPa to 3 GPa. Instead, I — Al and I + Al decrease and increase 
significantly, respectively. 

c. Reversibility. A first-order transition is frequently characterized by a two-phase co- 
existence region of metastable and stable phases as a function of the thermodynamic variable 
(e.g., temperature or pressure). In piezoelectric materials this is one source of irreversibil- 
ity (other significant contribution being due to the irreversible domain wall motion). It 
is interesting to note that the recovery run, carried out after the high-pressure and high- 
temperature cycles, revealed that the rhombohedral phase fraction had increased when com- 
pared to the prior the high-pressure situation. This suggests that high-pressure synthesis 
is a useful way to prepare single-phase rhombohedral ceramics in the vicinity of the MPB. 
The advantage over the Zr-rich rhombohedral ceramics is that in the vicinity of the phase 
transition certain piezoelectric constants are more susceptible to external stimuli. We note 
that recent neutron powder— and single crystal^ diffraction studies revealed that there is a 
secondary monoclinic Cm phase present in the Zr-rich case, together with the rhombohedral 
R3m/R3c phases. Recent single crystal study also showed that the diffraction data, collected 
on Pb(Zr .54Ti .46)O3 and Pb(Zr .69Ti . 31)03 samples are better interpreted in terms of the 
rhombohedral and monoclinic phases, rather than by the adaptive phase model^. The two- 
phase co-existence and the nature of the phase transition are believed to be crucial for the 
piezoelectric properties. 

IV. CONCLUSIONS 

High-pressure neutron powder diffraction experiments were applied to the classical piezo- 
electric compound, Pb(Zr .54Tio.46)0 3 . Weight fraction changes between the rhombohe- 
dral R3c and monoclinic Cm (low-pressures and room temperature) or between tetragonal 
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FIG. 6. (a) 5-cation (Zr or Ti) and oxygen bond lengths in the rhombohedral phase. The difference 
between B — Oai+i and B — Oai-i bond lengths increases with increasing pressure. The decrease 
in difference seen at 1.41 GPa pressure is probably related to the vicity of the transition to the 
cubic phase, (b) The distance D(l + Al, I — Al) between the oxygen triangles. In both panels, the 
3 GPa data are indicated by dotted lines. The inset shows the displacement of the B cations under 
pressure. At ambient pressures the B is closer to the larger triangle and displaces towards smaller 
triangle under pressure. 



PAmm phases as a function of hydrostatic pressure and function were determined. The 
Cm phase was observed only at low-pressures and ambient temperatures as it transformed 
to the PAmm phase at approximately 1 GPa and 400 K. As the earlier computations pre- 
dicted, the rhombohedral phase was favored at higher pressures, whereas the added heat 
increased the monoclinic phase fraction at constant pressure. This largely contributes to the 

12 



extrinsic piezoelectricity. These findings are in line with the computational model according 
to which the phase boundary between the rhombohedral and tetragonal phase in pressure- 
temperature plane is dictated by the two competing terms, octahedral tilting and entropy 
term. No support for a continuous polarization rotation was found. The oxygen octahedra 
was significantly distorted under pressure, accompanied by a significant displacement of the 
B cations. This contributes to the intrinsic piezoelectricity. After the experiments the frac- 
tion of the R3c phase was larger than initially, suggesting that the Cm phase is not stable. 
This is consistent with the first-order phase transition Cm — > R3c. 
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